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Abstract 

Given a foliation S of a manifold M , a distribution Z in M transveral 
to S and a Poisson bivector It on M we present a geometric method of 
reducing this operator on the foliation <S along the distribution Z. It 
encompasses the classical ideas of Dirac (Dirac reduction) and more mod- 
ern theory of J. Marsden and T. Ratiu, but our method leads to formulas 
that allow for an explicit calculation of the reduced Poisson bracket. More- 
over, we analyse the reduction of Hamiltonian systems corresponding to 
the bivector II. 
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1 Introduction 

The reduction theory of dynamical systems consist of two branches: the first 
branch deals with constrained Lagrangian systems, the second with constrained 
Hamiltonian systems. In the Lagrangian approach one consider separately the 
case of holonomic constraints, i.e. the constraints which may depend on veloci- 
ties, but only in such a way that the equations of constraints can be integrated 
to eliminate velocities, and the non-holonomic case. In many of these papers 
one first considers the Lagrangian formulation and then passess to the corre- 
sponding Hamitlonian formulation (see for examlple jj). The reduction theory 
in the Hamiltonian context has been initiated by P.A.M. Dirac, who in his fa- 
mous paper [3j described a method of reducing a given Poisson bracket onto a 
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submanifold given by some constraints <p provided they were of "second class". 
In this approach the classical notion of holonomic constraints is usually not 
introduced as in this context there is no obvious division of variables between 
"position" and "velocity" (or "momenta"). Recently, there has been much in- 
terest to extend the theory onto the case of generalized Hamiltonian systems 
(see for example [3]). 

The ideas of Dirac were developed in many papers, among others in P]-|10j 
(see also the literature quoted there). A geometric meaning of this reduction 
procedure has been investigated in |S] and in [TT] . 

In this paper we develop the ideas of |2j and [S] and present a constructive, 
computable method of reducing (locally) a given Poisson operator II to any 
regular submanifold S. The idea of the method is to choose a distribution Z 
(not necessarily integrable) that is i) transversal to the foliation S ii) at any 
x G M it completes T X S to T X M. iii) it makes the operator II Z-invariant (see 
definitions below) and then to deform the Poisson operator II to a new Poisson 
operator such that its image will be tangent to the submanifold S. This 
new operator will be always Poisson (and so its natural restriction to S will 
be Poisson). In consequence, we obtain a method of reducing a Hamiltonian 
system on M to a Hamiltonian system on every leaf <S„ of the foliation <S. This 
reduced system strongly depends on the choice of the distribution Z. As a 
special case we obtain the classical Dirac reduction of Hamiltonian system. All 
our considerations will be local in the sense that our manifold Ai is perhaps 
only an open submanifold of a larger manifold. Our construction is equivalent 
to the reduction method proposed by Marsden and Ratiu in 8 . However, our 
approach has advantages: it can be performed simultaneously on any leaf S v of 
the foliation 5, it is constructive (the approach of Marsden and Ratiu requires 
calculations of prolongations of Z-invariant functions and as such is difficult to 
perform in practice) and it is formulated in the language of Poisson bivectors 
rather than Poisson brackets. 

We want to stress that this method does not require the submanifold to 
be given by a holonomic constraints on some configuration space. In fact, we 
do not require our manifold to be a cotangent bundle to any configurational 
manifold at all, but of course our construction covers this special case as well. 
For example, it covers the case discussed in pp, where the auhtors obtain the 
Poisson operator on the constrained submanifold only in the case of holonomic 
constraints - because they simply restrict their Poisson operator n onto the 
constrained submanifold, and such restriction usually (apart form the holonomic 
case) destroys the poissonity of the operator. 

We have to notice that some particular versions of the proposed scheme 
appeared recently in ^2] an d Q3|- I n QH the author applied the same setting 
as above but with no conclusive formulas for computing the actual deformed 
Poisson bracket and since he did not use the notion of Z-invariance, his 
deformed operator Hd, although formally identical with our construction, was 
not always Poisson (it has been called there a pseudo-Poisson operator). 

Some basic steps of our construction have been presented in our previous 
paper ^Hj- This paper, however, was mainly devoted to completion of the 
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above picture by its "dual" part by developing a theory of a Marsden-Ratiu 
type reduction of presymplectic 2-forms SI that are (in a sense) dual to a given 
Poisson operator II. The picture presented here is much clearer and moreover 
it is paramietrization-free in the sense that we prove the main results without 
neccessity of dicsussing some particular functions defining our foliation S. More- 
over, in the paper JS] we focused on Dirac case while this paper has a general 
character. In the end, in this paper we also consider the related reduction of 
Hamiltonian dynamics. 

2 Geometric reduction of Poisson bivectors 

Let us consider a smooth manifold M. of arbitrary (finite) dimension n and a 
foliation S of M. consisting of the leaves S u parametrized by v G H k (so that 
k G N is the codimension of every leaf <S„). Consider also a regular distribution 
Z on M. (that is a smooth collection of the spaces Z x C T X M. where v is such 
that x G S v ) such that it completes every T X S to T X M in the sense that 

T X M =T X S V © Z x (1) 

for every x inAL Here and in what follows © denotes the direct sum of vector 
spaces. It means that every vector field X on the manifold M. has a unique 
decomposition X = Xy + X± such that for every x in M. the vector (X||) G 
T X S (X\\ is tangent to the leaves of the foliation S) while (X±) x G Z x (X± is 
contained in the distribution Z). The splitting induces the following splitting 
of the corresponding dual space T*A4: 

T*M =T* X S U © Z* (2) 

where T*S U is the annihilator of Z x while the space Z* is the annihilator of 
T X S V . Thus, any one- form a on M. has a unique decomposition a = an + aj_ 
such that (ay) G T*S (ay annihilates the vectors from Z) while {ot±) x G Z* 
(a± annihilates the vectors tangent to the foliation S ). We will call Xn and an 
as projections of X and a (respectively) on the foliation S. Abusing notation a 
bit we will write that X C TS if X = X\\ , X C Z if X = X± and similarly for 
one forms: a C T*S if a = an, a C Z* if a = 

Let us now suppose that our manifold M. is equipped with a Poisson bivector 
II (i.e. a bivector with vanishing Schouten bracket, sec 6 ). This operator 
induces the Poisson bracket {F, G} n = (dF, HdG) on the algebra of smooth 
functions on M., where (•, •) is the dual map between T*M. and TM.. A smooth 
real- valued function F on M. is called Z-invariant if the Lie derivative LzF = 
for any vector field Z C Z. We will now adopt the following definition. 

Definition 1 The operator IT is said to be Z-invariant if Lz {F, G} n = for 

any pair of Z-invariant functions F and G and every vector field Z C Z. 

Notice that our definition does not necessarily mean that LzII = for all 
vector fields Z C Z, as for any pair F, G of Z-invariant functions the condition 
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Lz {F, G} n = means only that the function (dF, (L^II) dG) vanishes. Thus, 
II does not have to be an invariant of the distribution Z to be iJ-invariant in our 
meaning. Notice also, that the above definition is equivalent to the statement 
that for any pair a, ft C T*S we have (a, (L^II) P) = (since if F is Z-invariant 
then dF C T*S). 

Suppose for the moment that the distribution Z is spanned by k vector fields 
Zi. We say, that the operator II is Vaisman |16| with respect to Z if for every 
vector field Zi there exists vector fields Wij, j = 1, . . .k such that 

L Zi n=J2W ij AZ j . (3) 

3=1 

It is easy to see that this definition does not depend on the choice of basis in 
Z (although the vector fields Wij obviously do). If the operator II is Vaisman 
with respect to Z, then it is also Z- invariant, as then for any two one- forms 
a,p C T*S 

(a, (L Z U) p) = £ (a, (Wij A Zj) p) = 0, 

since a and /3 annihilate all the vector fields Zi. The converse statement is 
however not true in general. 

Let us now consider a Poisson operator II on M. and define the following 
bivector: 

IId (a, p) = n(aii, P\\) for any pair a, P of one-forms. (4) 

We will often call the bivector II a deformation of II. This bivector occurred 
for example in I12| in a more restrictive context. Observe that it always exists 
and that it is uniquely defined once the foliation S and the distribution Z are 
given (it is thus a purely geometric construction) . It has an important property: 
its image lies always in TS. 

Lemma 2 Hjy(a) C TS for any one-form a on M., i.e. the image of Hp is 
tangent to the foliation S. 

Proof. We have to show that (P, UdCx) = for any p C Z*. But 

(P,H D a) =U D (p,a) =H(p h a\\) =0 

since /3|| = for every /3 C Z* . m 

Thus, the deformed bivector Ud has its image in TS and if we consider it as 
mapping from one-forms to vector fields on A4 then it can be naturally restricted 
to a bivector -k^ v on every leaf S v of S by simply restricting its domain to S v : 

n R« = ^ D \s v ■ 
Moreover, it induces a new bracket for functions on M.: 

{F,G} U =U D (dF,dG)=U((dF), l ,(dG),X 
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Of course, the bivector Up (and thus even ttr v ) in general does not have to be 
Poisson. However, it turns out that if II is Z- invariant then Hp (and thus every 
7Tfl„) is Poisson. 

Theorem 3 If II is Z-invariant then Hp given by ^ is Poisson. 

Proof. We will prove that the bracket {•, -} n is Poisson. Obviously, this 
bracket is antisymmetric and satisfies the Leibniz property. It remains to show, 
that it also satisfies the Jacobi identity, that is {{F, G} Ud , -^} n + cvc l- = 0, 
for any functions F,G,H. Using the definition of Hd, this condition can be 
written as 

(d((dF)u ' n ( dG )||))n ,n(dH)„^+qrcL = 0. (5) 
However, for any vector field Z C Z we have 

(d((dF)||,n(dG)||),z) = z(((dF)|,,n(dGO||)) =L z ({dF) v U{dG)^=Q 

due to the assumed Z-invariance of n. This means that diAdF)^ , H(dG)^) C 
T*S, so that 

(d ( (dF) „ , H (dG) |, )) = d ( (dF) || , n (dG) „ ) , 

and thus condition J5j turns out to be the Jacobi identity for n, which is satisfied 
since n is Poisson. ■ 

Thus, given a foliation S on A4 and a transversal distribution ZonM (such 
that |QJ is satisfied) we can reduce any Poisson bivector n that is Z-invariant to 
a Poisson bivector ttr^ on the leaf S v of S by deforming II to lie and then by 
restricting IIb to S„. This construction yields the same operator itr^ as in the 
approach of Marsden and Ratiu [S] ■ We will however show how this construction 
can be easily realized in practice. 

Remark 4 In case when the foliation S coincides with the symplectic foliation 
o/n we have of course that Hd = n , since in this case Tl((dF)j_) = for any 
function F and so {F, G} Ud = U((dF) {{ , (dG) [{ ) = U (dF, dG) = {F,G} n . In 
this case ttr v is the standard reduction of H on its symplectic leafS u . 

Let us now consider some special cases of our general situation. We firstly 
observe that the annihilator Z* of TS is defined as soon as the foliation S is 
determined (we do not need to specify a particular Z in order to define Z*). 

Definition 5 The distribution T> = II (Z*) (so that T> x = LT(Z*) ) is called a 
Dirac distribution associated with the foliation S. 
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Thus, the distribution T> is determined by S and by II. A priori, two limit 
cases are possible here. If TM =V(B TS we say that we are in the Dirac case. 
If V C TS we say that we are in the tangent case. In the Dirac case we have 
a canonical choice of Z: we can choose Z — T> (in this case IT is automatically 
Z-invariant, since Z is spanned by the vector fields Hamiltonian with respect 
to II). Nevertheless, we can also choose some other distribution 2 / D. In 
the tangent case we have no canonical choice of Z and we are free to find a 
distribution Z that makes IT Z-invariant. Anyway, in both cases we have many 
non-equivalent deformations IT^ (and thus projections -kr„). Generically, the 
distribution V will not be tangent to S, but it will not suffice to span TA4 
together with TS. 

Let us now suppose that the foliation <S of M is parametrized by the set of k 
functionally independent real valued functions ipi(x) so that its leaves have the 
form S v = {x G M. : ipi{x) = f,, f, € R, i = 1, . . . , k } where k is - as above - 
the codimension of the foliation. We will show how the above considerations can 
be written in the parametrization {ifi}. The one-forms difi constitute a basis 
in Z*. Then, the Dirac distribution V is spanned by k (possibly dependent) 
Hamiltonian vector fields A, = Udipi. Let us denote a basis of Z dual to the 
basis {d(fi} in Z* by Z- L , so that Zi((pj) = Sij. Our projections A|| and ay are 
then given by 

k 

X\\=X-'£X{<p i )Z i , 

i=l 

(obviously X\\((pi) =0 for all i so that indeed this vector field is tangent to the 
leaves of S) and by 

k 

a \\ = a ~ E a (^)^> 
i=l 

(and obviously ot\\(Zi) = for all i). Thus 

I%l|,/?||) = n (a- £a(Zi)dtpi,l3- ZPiZ^chpj 

= II(a,/?)- Y,p{Z j )TUa,dtp i )- Ea(Z i )n(rf^,/3)+ 

j = l i=l 

+ £ a(Z i )/3(Z j )n(# i ,rf^), 
so that the deformation Ud can be expressed as 

ITz, = II - £ A, A + i X>y^ A Zj (6) 

i £ i,j 

where the functions tpij are defined as 

ifij = {(Pi,yj} n =Xj(ipi) =U(dipi,d(pj). 

In the Dirac case all the vector fields Xi are transversal to the foliation S and 
are moreover linearly independent. It happens precisely when det((fij) ^ (the 
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functions ipi are then 'second class constraints' in the terminology of Dirac).The 
vector fields Zi (the dual basis to {dipi}) can be expressed through the vector 
fields Xj as 

k 

i=i 

Indeed, 

k k 

z Avi) = E {v~ l )sjX s {ipi) = £ (v3~ 1 ) si ^ s = <%. 

s=l s=l 

Moreover, in this case the deformation © attains the form: 

k 

n D = n-i^x j A4 (7) 

i=l 

This operator defines the following bracket on C°°(A4) 

k 

{F, G}n D = {F, G} n - £ fft}n(v _1 )« G} n , (8) 

which is just the well known Dirac deformation |2] of the bracket {., .}n- 

In the tangent case all the vector fields X{ are tangent to the foliation S 
and the deformation JfjJ) attains the form 

fc 

U D = U-J2X i AZ i , (9) 
i=i 

and has been considered in ^21 an d in 

In we have considered a bit more general then © deformation of II 
of the form lid = II — J2i ^ ^ ^i- The image of 11^ does not have to lie in 
TS. However, one can prove that it happens precisely when the vector fields Vi 
satisfy the following functional equation 

Vi=Xi + Y, V M) Z r ( 10 ) 

j 

and in this case 11^ is Poisson. Substituting this formula into 11^ yields 

n d = n - Y,Xi AZi+ £ v^z, a z h (11) 

and it can be proved (after some technical manipulations) that in this case 
n<j = Ud (even though V} (</?,) does not have to be equal to \ipij). A natural 
solution of the equation (|10|) in the Dirac case is given by Vi = \Xi and in the 
tangent case by Vi = Xi, which turns into the deformations (7J and 10 
respectively. Thus, in a sense, the deformation Hp is the canonical deformation 
in our setting. 
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Our construction is strongly related with the construction of Marsden and 
Ratiu. J. Marsden and T. Ratiu presented in |S] a natural way of reducing 
of a given Poisson bracket {•, -} n on M to a Poisson bracket {■,•}„. on a 
given submanifold 5„ (in our notation). Their method is non-constructive in 
the sense that in order to find the bracket {/, g}^ of two functions f,g : S v — ► 
R one has to calculate Z\s„ -invariant prolongations of these functions. Our 
construction is performed on the level of bivectors rather than on the level of 
Poisson brackets. This construction (by deformation of the bivector IT) applies 
directly to every leaf of the distribution S and moreover it is constructive. At 
every leaf, however, both constructions are equivalent. On the other hand, we 
make the assumption about the transversality of the distribution Z that was 
not present in the original paper of Marsden and Ratiu. This assumption is 
however very natural since it makes all the assumptions of Poisson Reduction 
Theorem in [H] automatically satisfied. 



3 Reduction of Hamiltonian dynamics 

Let us begin this section by stating - in our setting - a well-known theorem about 
the relation between the Dirac deformation IId of II and the dynamic imposed 
by the constraints. Suppose, thus, that our manifold M is a cotangent bundle 
of a Riemannian manifold Q with a covariant metric tensor g, so that M. — 
T*Q. Denote the corresponding contravariant metric tensor by G. Considar a 
Lagrangian dynamical system on TQ 

d dL dL n 

— — = 0, i = l,...,n, 12 

dt dq . % 

.t . 

with a potential Lagrangian function L(q,q) — ^q gq — V(q). This of course 
leads to a Hamiltonian equations of motion on M. = T*Q 

q i = {q i ,H} n , p t = {p l ,H} u (13) 

with the Hamiltonian H — ^p l Gp + V(q) and with the canonical Poisson oper- 
ator n. Let us now impose a physical constraint ip(q) — on our system and 
assume that in the beginning the coordinates of the system lie on the subman- 
ifold Qo of Q given by ip(q) = 0. One often makes a physical assumption here 
that the surface Qo starts to act on our system with a reaction force R(q, q) that 
is orthogonal to Qo and such that the trajectorties of the constrained system 

d dL dL ■ 
dt dq . % 

that start on Q remain on Q . On the level of the phase space M = T*Q the 
Hamiltonian system l|13f) is now subordinate to a pair of constraints: 

pi(g) = <p(q) = 0, <p 2 (q) = (V^) 4 Gp = 0, 
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where \7ip is the gradient of <p with respect to g-variables (differential of ip in 
Q), (the second constraint is a consequence of the fact that the velocities q must 
remain tangent to Qq) and thus modifies to 

<Zi = fe^On , Pi = {Pi,H} n + Ri{q,p). (15) 

Theorem 6 The equations U5\) are Hamiltonian and can be written in the form 

Qi = Ui,H} UD , p. l = {p l ,H} nD (16) 

where Tip = II — i Y^ii=i ^-i A Zi is the Dirac deformation of H given by the 
constraints ipi and ip2 ■ 

Thus, the response of the Lagrangian system l|12|) subordinated to the reac- 
tion forces R can be accounted for by the corresponding Direac deformation of 
the Poisson operator II. We will only sketch the proof. 

Proof. The reaction force R can be calculated by differentiating the as- 
sumed identity tp{q(t)) = twice with respect to t and eliminating the second 
derivatives qi with the help of equations (jl4|) and by using the demand that the 
force should be orthogonal to Qq. After some calculations we obtain that: 

R(q, P) = ■= , t * ({VpfGVV- (p t G)H v (Gp) +A)Vp (17) 
(Vp) G v<p K ' 

ri 2 

where H v is the Hessian of ip: (H^)^ — d °.g and where A = A(q,p) is given 
by 

s t/ys r,m,i,j 

so that it vanishes in the Euclidean coordinates when all Christoffcrs symbols 
Tj k are equal to zero. On the other hand, calculating the explicit form of l|16fl 
on the submanifold of M. = T*Q given by the constraints <^i,<£>2 leads to the 
equations (fT5|l with R given by (fTTjl . ■ 

Let us now consider a Hamiltonian vector field X = HdH on a general Pois- 
son manifold M. where H is some real-valued smooth function on M. (Hamilto- 
nian function). We constantly assume that we have a smooth, regular foliation 
S on M. and a regular distribution Z on M. such that Q is satisfied. The 
corresponding IId defined by (@J is Poisson and has its image tangent to the 
foliation S, so that it can be properly restricted on every leaf cv of S. Thus, 
the following definition makes sense. 

Definition 7 We call the vector field X o = HdcIH the Hamiltonian projection 
of the Hamiltonian vector field X = WdH . 

The vector field Xu lives on every leaf of the foliation in the sense that its 
restriction on the leaf S v is tangent to S v . Moreover, on the leaf <5>„ it coincides 
with the Hamiltonian vector field iTR u dh: 

Ht)dH\ s = TTji v dh, 
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where h = H\ s is the restriction of the Hamiltonian H to the leaf S v . To 
see this it is enough to choose a parametrization {ipi} of S and to pass to any 
system of coordinates of the form (x, ip). In these coordinates the bivector LTd 
has a matrix form with a non-zero upper-left block coinciding with the matrix 
form of tvr v and with the remaining terms equal to zero. 

There is a connection between X = TidH and its Hamiltonian projection 
X D = Tl D dH. 

Theorem 8 (Dynamic reduction theorem) If X = TldH , Xd = TlodH and 
Xi = Hdifi then 

X D =X ]] -^2Z i (H)X it{] . (18) 
i=i 

Proof. A direct calculation yields 
X D = U D dH = X-*t (Zi(H)Xi - Xi(H)Zi) + £ tpyZjfflZi = 

i=l i,j=l 

= X|, - Y,Zi{H) [Xi - jzVjiZj 
i=l \ J=l 

the last equality due to Xi(H) — (dH,Ud<pi) — — (dipi,HdH) = —X(ipi). Since 
<Pji = Xi(Lpj) it yields JTHJ- ■ 

k 

Observe that the difference between Xr> and X\\ is the term Zi{H)X i ^ 

i=i 

that is tangent to the foliation S, as it should be. 

Since for the Dirac case Xo = (by definition) we have 

Corollary 9 In the Dirac case Xp = Xu , so that in this case the Hamiltonian 
projection is just the natural projection (in the sense of direct sum) along the 
distribution 2. 

The term Xu in Xd has a well-known physical interpretation: it describes 
the evolution of the system X = HdH imposed with the constraints given by 
<Pi. The physical meaning of the second term in Xd is not clear for us, although 
it should represent some additional force (friction) acting on the system and 
tangent to the constraints. The authors are grateful for any hints in this matter. 

We are now in position to discuss the degeneracy of IId using the above 
Dynamic reduction theorem (Theorem [SJ. Let us first discuss the Dirac case. 

Proposition 10 Consider the Dirac deformation IId given by of a Poisson 
operator II on M. Suppose that the real valued functions a, i = 1, . . . , s on M. 
are such that they span the kernel of the operator H (i.e. are Casimir functions 
o/n in the sense that Hdci = 0) and such that the functions {c,, ifj} constitute 
a functionally independent set. Then 

i) the constraints <pi and the 'old' Casimir s Cj are all Casimir s o/IId 

ii ) any Casimir of IId rnust be of the form C{c.\ , . . . , c s , (pi , . . . ip^) ■ 
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Proof. The proof of i) is just a calculation. To prove ii), let us complete the 
functionally independent functions {c%)<Pj} to a coordinate system {cj, tpj, Xk} 
on M. Suppose that a function C — C(c,tp,x) is a Casimir of I±£> given by 
(0), i.e. that U D dC = 0. Then, according to Theorem [HI (HdC),, = 0, i.e. 
IliiC C Z. In the Dirac case the distribution Z is spanned by the vector 
fields Xi so that there must exist functions a; such that HdC = Y^t=i a i^i = 

SLi^n^i = ^-(l2i=i a i ( ^Pij- Thus, II ^dC- X)i=i a i^i) — or dC — 

J2i=i a id ( Pi + Z)i=iArfc» which proves ii). ■ 

Thus, we can state that the Dirac deformation (JJJ preserves all the old 
Casimir functions and introduces new Casimirs ifi and that no other Casimirs 
arise in this process. The situation in the general case is more complicated, 
since the Casimirs of II does not have to survive the general deformation J3J 
(or ©) and moreover since new Casimirs, different from tp^ ones, can arise. We 
can merely state that in the general case the function C is a Casimir of if 
and only if the vector field Y = HdC satisfies the relation 

i=l 

which for the Dirac case degenerates to the already discussed condition Yji = 0. 

4 Example 

Let us conclude this article by an example. Consider the so called first New- 
ton representation of the seventh-order stationary flow of the KdV hierarchy 
[TB ] .[T5 ] .[T7 j . It is the following system of second order Newton equations 

91, tt = -IO91 2 + 4<j 2 

92, m = -I69192 + IO91 3 + 493 

93, tt = -2O9193 - 8q 2 2 + 3O9192 - 15qf 

(where the subscript , t denotes the differentiation with respect to the evolution 
parameter t). By putting p\ — q^t, Vi — 92,tj P3 = 9i,t it can be written in a 
Hamiltonian form: 

d t 
^7(91,92,93,^1,^2,^3) =X = HdH 

where II is the canonical Poisson operator on the space M. = { (91 , 92 , 93 , Pi , P2 , Pz ) } , 



n = 






h ' 


. - J 3 






and with the Hamiltonian 

H = pip 3 + \p\ + 109i 2 93 ~ 4q 2 93 + 89192 2 - 109i 3 92 + 39i 5 . 
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Consider also a foliation S given by a pair of constraints 

</?l = ?3 + 9l<72, </?2 = Pi + qiP2 + 92P3- 

where ip2 is the so called G-consequence of tp\ (i.e. a lift from the configuration 
space {(qi, q2, 93)} to M) with respect to the antidiagonal metric tensor G |17| . 
The vector fields X; = Ildipi have the form 

9 9 d d d d d d 

Xi = -q 2 j. giT; -k— , X 2 = — + q 1 - Vq 2 j. p 2 j. P3^~ 

opi 0P2 0P3 oqi oq 2 oq 3 dpi dp 2 

and they are transversal to S so that we have the Dirac case. Thus, the distri- 
bution Z = T> — Sp{Xi} makes IT ^-invariant. The basis in Z that is dual to 
{dxpi} is Zi = ^X 2 , Z 2 = - ^Xt, where tp 12 = {ipi,<p 2 } u = 2q 2 + q\. The 
Dirac deformation ITd given by 10 attains in the adapted coordinate system 
{(qi, 32, <Pu<P2>P2iPa)} tnc form 

-9i -1 " 
2q 2 -q x 



P3 

-Pa 

It has, as it should, two Casimirs cpi, <p 2 . We can now easily restrict lio to the 
operator wr v on S v . If we parametrize S v with the coordinates {{qi,q 2l p 2l p 3 )} 
(the constraints ipi, (p 2 are constant on every S u ) then 

-qi -1 

2q 2 -qi 
qi -2q 2 p 3 

1 qi -P3 

which, in accordance with the theory, is non-degenerate. Observe that this 
expression actually does not depend on the choice of the leaf S v in the foliation 
S. The Hamiltonian projection Xd = liodH attains on every leaf <S„ the form 

1 ( d d a 9 9 \ 

WRv dhji = ■ — 2 Oil- \-Oi 2 - ^P 2 - ^P3-^— 

2<72 + q{ \ oqi dq 2 dp 2 op 3 J 

where the functions on, f3i are some rather complicated polynomial functions of 
coordinates and the parameters v± of the leaf S v . 

Let us now choose another distribution iJ, for which II is also ^-invariant. 
Since the operator II has a very simple form, any pair of constant fields will span 
a distribution Z that makes II -invariant, since then the Vaisman condition J2J 

is trivially satisfied. Thus, let us take Z = Sp + ^7 ^ + a^} (observe 

that this distribution is integrable). We have now to change the basis in Z to 



n 



D = 



1 

2<Z2 + q\ 



















































91 


-2<?2 








1 


qi 









2?2 + qf 
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a new basis {Zi, Z^\ such that the condition Zi(ipj) = Sij is satisfied. A simple 
calculation yields 



7 9 



1 - 92 



-'12 



_d_ 

dpi 



d 

dp 3 



7-2 



1 



d 



1 - 92 V dpi 



d 

dp 3 



Now, the general deformation © defined by the above distribution attains in 
the coordinates {(91, q%, fx, ^ajPaiPa)} the form 



n 



D = 



92-1 



1-92 

1 91 

















92-1 






qi - p 3 



-1 

-9i 



P3 - 91 




and thus the restricted operator ttr^ on the leaf S u parametrized with the co- 
ordinates {(91,92,^2,^3)} is 



92-1 






1-^2 

1 9i 





92-1 


91 - pa 



-1 
-9i 

P3 - 91 





which is again non-degenerate. Again, this expression does not depend on the 
choice of the leaf <S„ in the foliation S. 

In the end, consider yet another distribution that makes II ^-invariant, 
namely Z — Sp j + , j . The appropriate basis of Z is given by 

d d d „ d 

dpi ' 



Zi 



9 7 
92 ^ h , Z 2 

dpi dp 3 



so that we have Zi(ipj) 



The general deformation © defined by the above 



distribution yields in the coordinates {(91, 92, tpi, <P2,P2,Pz)} 



n 



D = 



92-1 

































1 












































-1 











9i 














-9i 






so that the restricted operator 7Tr u on the leaf S v (again parametrized with the 
coordinates {(91,92,^2,^3)}) is 



92-I 




1 


-9i 






9i 
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and is degenerate this time. 

Thus, given a foliation S, by choosing different distributions Z we can obtain 
several different Hamiltonian projections of our original Hamiltonian system, not 
only just the well known Dirac reduction. 

5 Conclusions 

In this article we formulated a comprehensive, geometric picture of what is know 
as Dirac reduction and Marsden-Ratiu reduction of a Poisson operator II on a 
foliation S not related with the symplectic foliation of IT. As a consequence, we 
obtain a geometric method of reducing of any Hamiltonian system to a Hamil- 
tonian system on the foliation S. Any such reduction depends merely on the 
choice of the distribution Z along which the reduction take place. Thus, we can 
reduce Hamiltonian systems on M. to Hamiltonian systems on S in many non- 
equivalent ways. However, the procedure of finding appropriate distributions Z 
(i.e. those that make n Z-invariant) is non-trivial and non-algorithmic. 

References 

[1] A. J. van der Schaft, B. M. Maschke, "On the Hamiltonian formulation 
of nonholonomic mechanical systems", Rep. Math. Phys. 34 (1994), no. 2, 
225-233. 

[2] P. A. M. Dirac, "Generalized Hamiltonian Dynamics", Can. J. Math. 2 
(1950) 129-148. 

[3] G. Blankcnstein and A. J. van der Schaft, Symmetry and reduction in im- 
plicit generalized Hamiltonian systems, Rep. Math. Phys. 47 (2001) 57-100 

[4] J. Sniatycki, Dirac brackets in geometric dynamics, Ann. Inst. H. Poincarc 
Sect. A (N.S.) 20 (1974) 365-372 

[5] A. Lichnerowicz, Variete symplectique et dynamique associee a une 
sous-variete, C.R.Acad.Sci.Paris Ser. A-B 280 (1975), A523-A527 

[6] A. Lichnerowicz, Les varietes de Poisson et leurs algebres de Lie associees, 
J. Diff. Geom. 12 (1977), 253-300. 

[7] J. E. Marsden and A. Weinstein, Reduction of symplectic manifolds with 
symmetry, Rep. Math. Phys. 5 (1974) 121-130 

[8] Marsden and T. Ratiu, Reduction of Poisson manifolds, Lett. Math. Phys. 
11 (1986) 161-169 

[9] J. Grabowski, G. Landi, G. Marmo and G. Vilasi, Generalized Reduction 
Procedure: Symplectic and Poisson Formalism Fortschr. Phys. 42 (1994) 
393-427 



14 



[10] M. Flato, A. Lichnerowicz and D. Sternheimer, Deformations of Poisson 
brackets, Dirac brackets and applications, J. Math. Phys. 17 (1976) 1754- 
1762 

[11] K. Marciniak, M. Blaszak, "Dirac reduction revisited", J. Nonlin. Math. 
Phys. 10 (2003) 451-463, or arXiv.org/nlin.SI/0303014. 

[12] L. Degiovanni, G. Magnano, " Tri-Hamiltonian vector fields, spectral curves 
and separation coordinates", Rev. Math. Phys. 14 (2002) 115-1163. 

[13] G. Falqui, M. Pcdroni, "Separation of variables for Bi-Hamiltonian sys- 
tems", Math.Phys.Anal.Geom. 6 (2003) 139-179. 

[14] Ch.-M. Marie, "Various approaches to conservative and nonconservative 
nonholonomic systems", Rep. Math. Phys. 42 (1998) 211-229. 

[15] M. Blaszak, K. Marciniak, "Dirac reduction of dual Poisson-presymplectic 
pairs", J. Phys. A: Math. Gen. 37 (2004) 5173-5187. 

[16] I. Vaisman, "Lectures on the Geometry of Poisson Manifolds" . Progress in 
Math., Birkhauser 1994. 

[17] M. Blaszak, K. Marciniak, "Separability preserving Dirac reductions of 
Poisson pencils on Riemannian manifolds", J. Phys. A: Math. Gen. 36 
(2003), 1337-1356. 

[18] S. Rauch-Wojciechowski, K. Marciniak, M. Blaszak, "Two Newton decom- 
positions of stationary flows of KdV and Harry Dym hierarchies" Physica 
A 233 (1996) pp.307-330 

[19] Blaszak M, "On separability of bi-Hamiltonian chain with degenerated 
Poisson structures", J. Math. Phys. 39, 3213 (1998) 



15 



